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We present the analysis of X-ray spectral variability made on a sample of 7 Seyfert 1 bright galaxies, using XMM-Newton 
data. From the "XMM-Newton Science Archive" we selected those bright Seyfert 1 showing one or more prominent flares 
in their 2-10 keV light curves. For each of them we extracted spectra in 3 different time intervals: before, during and after 
the flare. We fitted them with a simple power law and then shifted a narrow Gaussian emission and absorption line template 
across the 2.5-10 keV data, in order to investigate the presence of line-like features with a confidence level greater than 
99%. Some highly significant features were detected in 3 out of 7 sources studied. In particular, the 3 sources, namely 
PG 1211-hl43, NGC 4051 and NGC 3783, showed the presence of a variable emission feature in the 4.5-5.8 keV band, 
characterized by an increase of its intensity after the flare peak. Because of the observed variability pattern, this feature 
seems to be ascribable to a reverbered redshifted relativistic component of the Fe K line. 
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1 Introduction 

Radio-quiet Active Galactic Nuclei (AGN), Seyfert galax- 
ies in particular, show extreme X-ray variability phenom- 
ena (Mushotzky, Done & Pounds 1993). These are charac- 
terized by large amplitude flux variations (AFx/Fx^l) 
on small timescales (At < 1 day). Such kind of phenomena 
are usually indicated with the term of "flare" although a rig- 
orous definition of flare doesn't yet exist in literature. The 
physical mechanisms which produce an AGN X-ray flare 
are not understood. In particular it is not clear where the 
flares do originate (e.g. Lightman & Rybicki 1979, Payne 
1980, Poutanen & Fabian 1999, Malzac & Jourdain 2000). 
Despite the present uncertainties on the flares trigger mech- 
anism, such phenomena are useful tools to investigate the 
innermost regions of AGNs (geometry, chemical composi- 
tion, black hole mass, etc.). In the context of a disk-corona 
model (e.g. Haardt & Maraschi 1991,1993) we expect in- 
tense X-ray continuum flux variations to be followed by a 
disk reflection component response. In particular, if the flare 
originates close to the central black hole it might be possi- 
ble to observe a reverbered relativistic Fe K line response 
(e.g. Fabian et al. 2000). The XMM-Newton MCG -6-30-15 
observation is the first one in which a reverbered relativis- 
tic Fe K line has been detected after a huge hard X-ray flux 
flare (Ponti et al. 2004). 

In the following we present a systematic study of X-ray 
spectral variability in response to a X-ray flare in a sample 
of bright Seyfert 1 galaxies observed with XMM-Newton. 



Table 1 Sources of the sample: column (1) duration of 
the flare; column (2) 2-10 keV count rate variation, with 
respect to the average count rate, between the maximum and 
the minimum of the flare. 



Name 


Flare 


At 


ACR(2-10 keV) 






(ks) 








(1) 


(2) 


NGC 3783 


I 


52 


0.65 




II 


80 


0.81 


PG 1211+143 


I 


12 


1.00 


TONS 180 


I 


6 


0.90 




II 


9 


0.60 


NGC 4051 


I 


6 


1.58 




II 


8 


0.88 


IZWl 


I 


6 


1.00 


PG 1448+273 


I 


4 


0.61 


AKN 564 


I 


2 


1.09 



e-mail: demarco@sissa.it 



The aim is to detect transient emission/absorption features 
in the Fe K band, produced in response to the flare emission. 

2 Selection of the sample 

The initial sample of sources was obtained cross-correlating 
the "XMM-Newton Master Log & Public Archive" and the 
"Veron Quasars and AGNs" catalog (Veron-Cetty, Veron & 
Gonial ves 2001) in order to select only Seyfert 1, 1.2 and 
1.5 (public data up to 2004 December). The resulting sam- 
ple was cross-correlated with the "ROSAT All Sky Survey: 
Bright Sources" catalog (Voges et al. 1999) in order to as- 
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PG lSll + 143 



WGC 4051 




Time (s) 

Fig. 1 PG 121 1+143 2-10 keV light curve with time in- 
tervals chosen for spectral analysis. 
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Fig. 2 NGC 4051 0.3-2 and 2-10 keV light curves with 
time intervals chosen for spectral analysis. 



NGC 37B3 



sociate to each source a mean count rate. We excluded ob- 
jects with ROSAT count rate < 0.40 c/s (which, assuming 
a standard spectrum with F = 1.7 and galactic absorption, 
corresponds to a mean 0.1-2.5 keV flux of ~ 1.3 x 10~^^ 
erg cm~^ s~^); this should avoid the inclusion in the final 
sample of too weak objects for which a time resolved spec- 
tral analysis would be impossible. Also observations lasting 
less than 10 ks were excluded in order to increase the prob- 
ability to study the emission properties before, during and 
after the flares. The final sample contains 67 observations. 
The 2-10 keV light curves were extracted from the EPIC pn 
camera data. We looked for those ones in which at least one 
flare event was present. The flare was defined as a count rate 
variation between the minimum phase and the peak phase > 
50% of the mean count rate registered during the entire ob- 
servation, i.e. (CRmax — CRmin)/(CR) > 0.5. The final 
sample consists of 7 sources and 10 flares. They are listed 
in Tab[ll some examples of flares are reported in Fig[l}|3l 
The MCG -6-30-15 and MKN 766 observations as wefl as 
one observation of NGC 405 1 were excluded because they 
are characterized by the superposition of a large number of 
flares, making hard to identify each of them. 

2.1 Spectral analysis 

The spectra of the 7 sources were extracted in three different 
time intervals: before, during and after each flare (hereafter 
indicated with the letters P, F and D respectively). Because 
of the limited duration of the observation in the case of NGC 
3783 we chose the following time intervals for spectral anal- 
ysis (see Fig. O: one time interval characterized by rela- 
tively low flux (PI), two time intervals during the first flare 
(Fl and Dl) and three time intervals during the second flare 
(P2, F2 and D2). Only EPIC pn data were used. The analysis 
was carried in the 2.5-10 keV band, because in this energy 
band we expect to observe spectral features associated to the 
Fe K line, and also to avoid warm absorber contamination. 
In the case of NGC 3783 we considered only the 3.5-10 
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Fig. 3 NGC 3783 0.3-2 and 2-10 keV light curves with 
time intervals chosen for spectral analysis. 



keV band because in this source the warm absorber compo- 
nent extends to higher energies. Each spectrum was initially 
modeled with a simple power-law. Then a narrow Gaussian 
emission/absorption line template was added to the model 
and shifted along the entire energy band at steps of 50 eV. 
At each step the fit improvement with respect to the simple 
initial model was computed in terms of variations. An 
emission/absorption Ax^ vs gaussian centroid-energy plot 
was produced in order to search for line-like features with a 
confidence level greater than 99% (i.e. Ax^ > 9.21). As an 
example the 2.5-10 keV spectra of PG 121 1+143 are shown 
in Fig.|4|and the corresponding Ax^ plots are shown in Fig. 
13 Once significant features were detected, Gaussian lines 
have been added to the model fixing the width of each Gaus- 
sian model to its best fit value. We looked for line-intensity 
and line-energy variations comparing the Gaussian normal- 
ization vs energy parameters contour plots in the different 
time intervals. Each contour plot was obtained by fixing all 
the added Gaussians to their best fit values apart from the 
one for which we wanted to analyze the time variations (see 
Fig.|6Hi. 
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Fig. 7 NGC 405 1 : comparison among the 90% and 99% 
normalization vs. gaussian centroid energy contour plots for 
the variable feature at ~4.5-5.8 keV in the chosen time in- 
tervals. Fe K line intensity during periods Fl and Dl was 
significantly greater than during periods P2, F2 and D2. 



Fig. 4 PG 121 1+143: the 2.5-7.5 keV spectra during time 
intervals P, F and D. The EPIC pn data are used. 
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Fig. 5 PG 1211+143: Ax^ vs. emission gaussian 
centroid-energy plot in the chosen time intervals. 
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Fig. 6 PG 1211+143: comparison among the 90% and 
99% normalization vs. gaussian centroid energy contour 
plots for the variable feature at ^5.0-5.5 keV in the cho- 
sen time intervals. Fe K line intensity during period D was 
significantly greater than during periods P and F. 



^ NGC 37S3 




5 5.3 5.4 

LineE keV 

"tr, NGC 3783 



^ -Flare II 

t ; P2 




5 5.2 5.4 

LineE keV 



Fig. 8 NGC 3783: comparison among the 90% and 99% 
normalization vs. gaussian centroid energy contour plots for 
the variable feature at ~5. 1-5.5 keV in the chosen time in- 
tervals. Fe K line intensity during periods Dl and D2 was 
significantly greater than during periods PI, Fl, P2 and F2. 



3 Results and discussion 

Only 3 (PG 1211+143, NGC 4051 and NGC 3783) out of 
the 7 studied sources show highly significant emission and 
absorption features. Their light curves are showed in Fig[I]- 
[3l Moreover, these 3 sources display variable emission fea- 
tures in the 4.5-5.8 keV band, the intensity of this compo- 
nent rising after the flare peak (see Figl^HHl). This suggests 
a delayed response with respect to the continuum emission. 
The variability of the feature is detected with a confidence 
level greater than 99%. 

The hypothesis that the feature might be a curvature of the 
spectrum produced by variations of the warm absorber ion- 
ization parameter (^ = Lx/nR^, where n is the density of 
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the absorbing gas and R is its distance from the source) 
seems to be ruled out by the observed variability pattern: 
in fact, in this case we should expect a decrease of the spec- 
trum curvature intensity (increase of ^ and, hence, decrease 
of absorption) after the flare peak. Moreover, the fact that 
we don't observe contemporary variations of the narrow Fe 
Ka line core component makes implausible the possibil- 
ity that the feature is due to a reflection continuum effect. 
Hence it is most likely ascribable to a reverbered redshifted 
relativistic component of the Fe ka line, produced in the 
inner regions of the accretion disk. 

4 Conclusion 

Time-resolved spectral analysis of 7 bright Seyfert 1 has 
been presented. These sources have been selected for the 
presence, in their 2.5-10 keV light curve, of at least one 
prominent flare. For each source and for each flare, we an- 
alyzed the data before, during and after the X-ray emission 
peak. A variable emission feature was observed in 3 up to 7 
sources spectra at energies of 4.5-5.8 keV. It shows an in- 
crease of intensity after the flare peak at the 99% confidence 
level. In our view the feature can be interpreted as a broad 
highly redshifted component of the Fe Ka line produced in 
the innermost regions of the accretion disk and, hence, sub- 
jected to relativistic effects. The variability pattern suggests 
that it may be produced by the reverberation of the hard X- 
ray flux continuum. Assuming a disk-corona model, from 
the measured delays it is possible to give a maximum limit 
to the distance between the corona and the accretion disk. 
The obtained values are: 32 r^ for PG 121 1-^143, 800 r^ for 
NGC 4051, 32 r^ and 76 r^ from the first and the second 
flare of NGC 3783 respectively. 
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